Exact Performance of Concatenated Quantum Codes 
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When a logical qubit is protected using a quantum error-correcting code, the net effect of coding, 
decoherence (a physical channel acting on qubits in the codeword) and recovery can be represented 
exactly by an effective channel acting directly on the logical qubit. In this paper we describe a 
procedure for deriving the map between physical and effective channels that results from a given 
coding and recovery procedure. We show that the map for a concatenation of codes is given by 
the composition of the maps for the constituent codes. This perspective leads to an efficient means 
for calculating the exact performance of quantum codes with arbitrary levels of concatenation. 
We present explicit results for single-bit Pauli channels. For certain codes under the symmetric 
depolarizing channel, we use the coding maps to compute exact threshold error probabilities for 
achievability of perfect fidelity in the infinite concatenation limit. 

PACS numbers: 03.67.-a, 03.65.Yz, 03.67.Hk 



I. INTRODUCTION 

The methods of quantum error correction Jl|, |^, [| have, 
in principle, provided a means for suppressing destructive 
decoherence in quantum computer memories and quan- 
tum communication channels. In practice, however, a 
finite-sized error-correcting code can only protect against 
a subset of possible errors; one expects that protected in- 
formation will still degrade, albeit to a lesser degree. The 
problem of characterizing a quantum code's performance 
could thus be phrased as follows: what are the effective 
noise dynamics of the encoded information that result 
from the physical noise dynamics in the computing or 
communication device? 

One could address this question by direct simulation 
of the quantum dynamics and coding procedure. How- 
ever, for codes of non-trivial size, this approach rapidly 
becomes intractable. For example, in studies of fault- 
tolerance Jpj one often considers families of concatenated 
codes ||, ||]. An iV-qubit code concatenated with itself 
I times yields an iV -qubit code, providing better error 
resistance with increasing I. For even modest values of 
N and £, simulation of the resulting 2^ N •'-dimensional 
Hilbert space requires massive computational resources; 
using simulation to find the asymptotic performance as 
£ — > oo (as required for fault-tolerant applications) is 
simply not on option. 

Instead, a quantum code is often characterized by the 
set of discrete errors that it can perfectly correct || . For 
example, the Shor nine-bit code jlj was designed to per- 
fectly correct arbitrary decoherence acting on a single bit 
in the nine-bit register. Typical analyses of such codes 
implicitly assume that the physical dynamics can be de- 
scribed by single-bit errors occurring at some probabilis- 
tic rate; if this rate is small (e.g. 0(p) for p< 1), the 
probability that these errors will accumulate into a multi- 



bit uncorrectable error is also small (e.g. 0(p 2 )). This 
type of leading-order analysis is limited to a weak-noise 
regime, and to error models strongly resembling the er- 
rors against which the code protects. Outside of this 
regime, these approximation methods fail to accurately 
describe the evolution of the encoded information. 

In this work we take a different approach to charac- 
terizing error-correcting codes, which leads to a simple, 
exact analysis for arbitrary error models. As suggested 
above, a code transforms the physical dynamics of the 
device into the effective dynamics of the encoded infor- 
mation. In section II we derive this transformation for 
arbitrary noise, and present a compact method for its 
calculation. 

In the case of identical, uncorrelated noise on individ- 
ual qubits, this notion becomes particularly natural: en- 
coding a logical qubit in several physical qubits yields an 
evolution less noisy than if the logical qubit had been 
stored, unencoded, in a single physical qubit. Thus a 
code acts a map on the space of qubit dynamics, map- 
ping the dynamics of a single physical qubit to the dy- 
namics of the encoded logical qubit. In section III 
show how to calculate this map, and in section IV 



we 
we 



use these maps to dramatically simplify the calculation 
of effective dynamics for concatenated codes when the 
physical dynamics do not couple code blocks. 

In section |v|, we restrict our attention to uncorrelated 
single-bit Pauli errors, and in section VI we calculate 
the exact performance of several codes of inter est un- 
der these error models. Finally, in section VII we use 



the coding maps to calculate the performance of certain 
concatenated codes, and find the exact threshold error 
probability for perfect fidelity in the infinite concatena- 
tion limit. These thresholds serve as important figures 
of merit for concatenation schemes, and for the codes 
considered here we find that the traditional approximate 
method s und erestimate these thresholds by up to 44%. 
Section VIII concludes, suggesting potential future ap- 
plications for these techniques. 
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II. DESCRIBING CODE PERFORMANCE 
WITH EFFECTIVE CHANNELS 

In this section we first describe error-correcting codes 
using a language that will facilitate the subsequent de- 
velopment. We will then present our method for exactly 
describing the effective dynamics of the encoded infor- 
mation. Though for clarity we restrict our discussion 
to codes storing a single qubit (sometimes called k = 1 
codes) all of the presented methods generalize naturally 
to codes storing quantum information of arbitrary dimen- 
sion. 

As an important preliminary, it can be argued ^. 
that all physically possible transformations taking quan- 
tum states p on a Hilbert space Ti to states p' on a Hilbcrt 
space Ti' may be written in the following form: 

p^ P ' = Y,A jP A) with 5>t^- = l (1) 



where the Aj are linear operators from Ti to Ti' and 1 de- 
notes the identity operator on Ti. Such transformations 
are called quantum operations or channels, and are nec- 
essarily linear, trace-preserving, and completely positive. 
It is easy to see that the composition of quantum opera- 
tions is also a quantum operation. (One also sees defini- 
tions requiring only J2j AjAj < 1, corresponding to the 
weaker requirement that a quantum operation be trace 
non-increasing rather than trace-preserving. However, 
the requirement of trace-preservation is better suited to 
our purposes here. See jj| for a discussion of the distinc- 
tion.) 



A. The Error-Correction Process 

The error-correction process, consisting of encoding, 
noise, and decoding, is depicted in Figure y we consider 
each stage in turn. An W-qubit code C uses a register 
of N qubits to encode a single logical qubit a|0) + /3|1) 
by preparing the register in the state a|0) + /3|1), where 
|0) and |1) are orthogonal states in the 2 JV -dimcnsional 
Hilbert space of the register. The codespace (i.e., the 
space of initial register states) is spanned by these two 
states. In what follows it will be convenient to describe 
states by density matrices: let the logical qubit be given 
by p Q and the initial register state by p(0). Writing B = 
|0)(0| + |1)(1|, the encoding operation £ : po — > p(0) is 
given by 



P (0) =S[p ] =B Po B\ 



(2) 



As B^ B — |0)(0| + |1)(1| = 1, £ is a quantum operation. 

After the encoding, the register state evolves due to 
some noise dynamics. In the setting of a quantum com- 
puter memory, the dynamics are continuous in time; as- 
suming evolution for a time t, we have p(t) = A/"t[p(0)] 
with A/f a quantum operation depending continuously on 
t. (For master equation evolution p — C[p], we have 
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FIG. 1: The error-correction process: a logical single-qubit 
state po is encoded in an A-qubit register as p(0). A noise 
process transforms the register to state p(t), which is then 
decoded to yield the logical single-qubit state p/. 



Aft 



) We will often omit the subscript t and sim- 



ply write M . In the setting of a quantum communication 
channel, the noise process is usually given by the discrete 
application of a quantum operation Af; thus if the trans- 
mitted state is p(0), the received state is A/"[p(0)], which 
we write as p(t) for consistency. 

After the noise process, an attempt is made to recover 
the initial register state p(0) from the current register 
state p(t) by applying a quantum operation TZ, which 
may be written as 

K[p(t)} = J2A j p(t)A] with 5>tA, = l. (3) 



As the initial state p(0) is known to be in the codespace, 
it is clearly more beneficial to return the state p(t) to the 
codespace than to do otherwise: lacking any other infor- 
mation, one could at least prepare the completely mixed 
state in the codespace |(|0)(0| + |1)(1|), yielding an aver- 
age fidelity of \ , rather than leaving the register outside 
the codespace, yielding a fidelity of 0. We will therefore 
restrict our attention to error-correction processes TZ that 
take all register states back to the codespace. (I.e., we 
assume no leakage errors during recovery.) 

With the above assumption, the post-recovery state 
TZ[p{t)] has support entirely on the codespace; thus it 
can be described by its restriction to the codespace, the 
logical single-qubit state pf such that £[pf] — TZ[p(t)]. 
Call T> = £^ o TZ the decoding operation (shown to be a 
quantum operation in Lemma 1 of Appendix |b|) : 

V[p(t)} = B^TZ[ P {t)]B = ]T B^A jP {t)A\B. (4) 

j 

With this definition, pf — V[p(t)\. We will consider the 
logical state pf as the outcome of the error-correction 
process, and therefore may say that the code is given by 
its encoding and decoding operations, i.e. C — (£,T>). 

To build intuition for the decoding operation T>, we 
note that for most codes considered in the literature (and 
all of the specific codes considered later in this paper) the 
recovery procedure is given in a particular form. First, a 
syndrome measurement is made, projecting the register 
state onto one of 2 W_1 orthogonal two-dimensional sub- 
spaces; let the measurement be specified by projectors 
{Pj}. After the measurement (whose outcome is given by 
j, the index of the corresponding projector), the recov- 
ery operator Rj acts on the register, unitarily mapping 
the subspace projected by Pj back to the codespace. For 
such codes, the recovery superoperator is given by (H) 



3 



with Ah 



R P 



and lZ[p(t)] is the expected state that 



results from averaging over syndrome measurement out- 
comes. 

For codes of this form, let {\0j), |lj)} denote the or- 
thonormal basis for the syndrome space projected by 



Pj such that Rj\0j) 

RjPj 



= |0) and Rjll, 



1 



Then 



|0)(Q. 



and using the expression for 

V given in (|4|) yields 
T>W)\ = T. j &R } P jP <t)P}R\B 

= EidQXOil + liXiiDpWdOjXoi + 11^(11) 

(5) 

Thus Pj — T>[p(t)] is the sum of the single-qubit density 
matrices that result from restricting p{t) to each of the 
syndrome spaces, with basis {|0j),|lj)} determined by 
the recovery operator. 

As an example, consider the bitflip code ||, a three- 
qubit code that protects against single bitflip errors. The 
bitflip code's encoding transformation is given by 



|0) i ^ |0> = |000), |1) » |1> = |111). 



(6) 



After the action of some error dynamics, the syndrome 
measurement then projects the register state into one 
of four subspaces: the codespace itself, and the three 
subspaces that result from flipping the first, second, or 
third bit of states in the codespace. The corresponding 
recovery operator simply flips the appropriate bit back, 
attempting to reverse the error. Thus the basis specifying 
the decoding operation is given by 




|000), |1 

1 100), |li 

1010), |1 2 

1001), |1 3 




(7) 



We will use the bitflip code as an example throughout 
this work. 



arbitrary noise Af acting on ./V-qubit states. By calculat- 
ing such a compact description, we may easily find the 
effective evolution of an arbitrary initial state po without 
explicitly considering the physical noise dynamics. As 
we now show, Q may be written as a 4 x 4 matrix with 
a simple interpretation. (See § for a full discussion of 
qubit channels represented in this fashion.) 

For each Pauli matrix a £ {I,X,Y, Z}, let (a)o = 
tr(crpo)- The Pauli matrices form a basis for qubit density 
matrices, and so the initial logical qubit po may linearly 
parameterized by its expectation values (er)o as follows: 

Po = \(I)oI + \{X) Q X + \{Y) a Y + \{Z) Q Z. (9) 

(As the trace of a density matrix must be 1 we will 
always have (/) = 1, but it will be convenient to in- 
clude this term.) Similarly, the final logical qubit pf 
may be linearly parameterized by its expectation val- 
ues (a) f = tr(apf). Thus the effective channel Q 
may be written as the mapping from the expectation 
values (a)o of po to the expectation values (a) f of 
Pf . Writing p = ((I)o,(X) ,(Y)o,(Z)o) T and p f = 
((I) f , (X) f , (Y) f , (Z) /) T , the linearity of Q allows it to 
be written as the 4x4 matrix such that pf = Gpo- The 
fidelity of a pure logical qubit p through the effective 
channel is then given by tr(poP/) = \po Pf — jPq T Qpo- 
Thus to fully characterize the effective channel Q we need 
only find the entries of its 4 x 4 matrix representation. 
(More generally, if the code stored a d-dimensional state 
rather than the two-dimensional state of a qubit, the log- 
ical density matrices po and p / would be expanded in the 
basis of the identity matrix and the d 2 — 1 generators of 
SU(d), and Q would be represented as a d 2 x d 2 matrix.) 

To find these matrix elements, we consider the encod- 
ing and decoding processes in more detail. Letting E a 
denote \£, [aj, the encoding transformation £ acts on po 
(given by (pi)) to prepare the initial register state 

p(0) - (J) Ej + (X) E X + (Y)qEy + (Z) E Z . (10) 



B. Calculating the Effective Dynamics 

The transformation pa — > pf gives the effective dynam- 
ics of the encoded information resulting from the physi- 
cal dynamics Af. Let Q be the map giving these effective 
dynamics: Pf = G[po\- From the above discussion, the 
effective dynamics are simply the result of encoding, fol- 
lowed by noise, followed by decoding, i.e. 



Q = VoM oE. 



(8) 



As Q is the composition of quantum operations £, Af and 
T>, it is itself a quantum operation. We may therefore call 
G the effective channel describing the code C = (£,£>) 
and physical noise dynamics Af. 

Because the effective channel Q is only a map on single 
qubit states, it should have a compact description — in 
particular, a description much more compact than some 



Thus the encoding operation £ is completely character- 
ized by the E a operators, which are easily constructed 
from the codewords: 



Ex 

Ey 

E z 



= |(|0)<0| 



(|0)(1| 



f |T)<1|) 
+ |T)<0|) 
T|+*|I}(0|) 



IHlo) 

|(|0)(0|-|T)(T|) 



(11) 



As expected, p(0) is the state po on the codespace, and 
vanishes elsewhere. 

Now consider the decoding process, which yields the 
logical state pf. We may express the expectation values 
(a) f in terms of p(t), the register state prior to recovery, 
as follows: 



(a) f = tv(<7 Pf )=tv(aV[p(t)}) 
= tr (V . aBtA jP (t)A]B 



(12) 
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Exploiting the cyclic property of the trace and noting 
that BoB^ =£[<j]= 2E a , we have 

(a) f = tr(D a p(t)) where D a = 2^A]E <J A j . (13) 



Thus the decoding operation T> is completely character- 
ized by the D a operators. 

Substituting p(t) = Af[p(0)] into (|l3|), we have (a) f = 
ti(D a Af\p(0)]). Substituting in the expression for p(0) 
given by (ulm then yields 



(a) f = tr D a Af 



(14) 



Letting the matrix elements of Q be given by 

Gaa' = tx{D a Af\E a ,\) (15) 

for a,<r' S {I,X,Y,Z}, we have (<7>/ = £ CT , G CT(7 ' (</)o, 
i.e. p f = Gpo- 

To completely characterize the effective channel G, 
then, we need only compute these matrix elements. In 
fact, trace-preservation (i.e. (I) / = (/)o) requires Gn = 
1 and Gix = Giy = Giz = 0. Thus the effect on the 
logical information of the potentially complex dynamics 
of the 7V-qubit register space are characterized by the 
remaining twelve matrix elements of G ■ If Af is time- 
dependent, then the only observable effects of this time- 
dependence will appear in the time dependence of the 
Gaa' , and Gt gives the effective channel for correction per- 
formed at time t. Note that the dynamics Af need not 
be related to those against which the code was designed 
to protect. 

We have thus shown that the effective dynamics may 
be calculated by evaluating (|l^), which requires con- 
structing the E„ and D a operators. The E a operators 
are easily understood to be the operators which act as 
■id on the codespace and vanish elsewhere; to build in- 
tuition for the D a operators, consider codes whose re- 
covery is specified by syndrome measurement projectors 
\Pj} and recovery operators {Rj} as discussed in sec- 
tion [II A . For these codes, we have Aj 



RjPj, and so 

D (j = 2J2j Pj RjE a RjPj. This expression may be sim- 
plified by noting that E a maps the codespace to itself 
and vanishes elsewhere, and Rj unitarily maps the space 
projected by Pj to the codespace. Thus R^E a Rj unitar- 
ily maps the space projected by Pj to itself and vanishes 
elsewhere, i.e. P d R}jE„RjP 2 = R^EaRj- We therefore 
have 



2^4^,, 



(16) 



Using the expressions for E„ given in (^lj) and Rj = 
|0> (0,1 + 11) (1,|, we have 

Di = E,m(o 3 \ + \i 3 )(h\) 

Dy = J2j(-i\0j){lj\+i\lj){0j\) 
Dz = EiQ0s)(0j\-\Wi\)- 



(17) 



Thus we see that in this case D a is simply the sum of the 
operators a acting on each of the syndrome spaces, with 
Z-eigenstates |0j) and \Xj) determined by the recovery 
procedure. Note that D] is the identity operator on the 
entire register space. 



III. CODING AS A MAP ON CHANNELS 

One often considers noise models Af consisting of un- 
corrected noise on each of the N physical qubits. This 
type of model arises naturally in a communication set- 
ting, where the register qubits are sent over a noisy trans- 
mission line one at a time, and is also appropriate for var- 
ious physical implementations of a quantum computer. 
(By contrast, one can also consider error models in which 
correlated noise dominates Q.) For such models, we may 
write 



(18) 



where Af^ is a quantum operation on a single qubit. 

The goal of encoding a qubit is to suppress decoher- 
ence: multiple qubits are employed to yield an effective 
channel G, which should be less noisy than the channel 
resulting from storing information in a single physical 
qubit, namely Af^. A code can thus be seen as a map 
on channels, taking TV*- 1 -* to G- More precisely, for an N- 
qubit code C = (£,£>), define the corresponding coding 
map ft c by 

n c :AA (1) -^G =VoAf {1) ® N < 



e. 



(19) 



We now derive an expression for the coding map f2 

we de- 



ll B 



of an arbitrary code C = (£,V). In section 
scribed how G may be specified by its matrix elements 
Gaa' , given by (|l5|). Since Af^ is a single-qubit quantum 
operation, it may also be written as a 4 x 4 matrix such 
that if W (1) takes p to pt , then p 1 = Af^p . We seek an 
expression for the matrix elements of the effective chan- 
nel G in terms of the matrix elements of the physical 
channel Af^. 

Operators on N qubits may be written as sums of ten- 
sor products of N Pauli matrices; we may therefore write 
the E a and D a operators describing C = (£,£>) as 



Ea' = a {li}(h^) 



{I ,X,Y,Z} 



V N . 



(W) (2°) 
(21) 



{I ,X,Y,Z) 



E.g., for the bitflip code described in section [II A| by (|6|) 
and (j^), we may calculate the E a > and D a operators using 
( |TT| ) and (^); expanding the results in the basis of Pauli 
operators yields 



Ei = |( III + IZZ 
E x = f ( XXX - XYY 
E Y = I( -YYY + YXX 
E z = f ZZZ + ZII 



ZIZ + ZZI ) 
YXY - YYX ) 
XYX + XXY ) 

IZI + IIZ ) 



(22) 
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and 



concatenated code: 



Di 
D x 
D Y 
D z 



III 

XXX 
|( YYY 
|( -ZZZ 



YXX 
ZII 



XYX 
IZI 



XXY ) 
II Z ). 



(23) 



To find the matrix elements of the effective channel, 
we substitute ((l^), (|20|), and ( f2l| ) into the expression 
for these matrix elements given by ([l5]). Noting that 
7V[i^i ® . . . <g> i/xjv] = A/" (1) ® . . . <8> A/" (1) and 
tr((4 <g> B)(C <g> £>)) = tr(AC)tr(BD) yields 

, w \ 
^}«U}II tr ^W[lMi]) ■ (24) 



From the orthogonality of Pauli matrices, the matrix 
when written as a vector of expectation values, has a 1 in 
the /ii component and zeros elsewhere. Further, tr(;/jp) 
is simply the V\ component of p. Thus tr^AA 1 ) [\p-i]) = 

i and we have 




Goo' 




(25) 



Thus the matrix elements of Q can be expressed as 
polynomials of the matrix elements of Af^ x \ with the 
polynomial coefficients depending only on the E a i and 
D a of the code. These polynomials specify fi ! . By com- 
puting these polynomials for a code C, one can easily 
calculate the effective channel for the code C due to any 
error model with identical, uncorrelated noise acting on 
each physical qubit. (If a different noise model acts on 
each physical qubit, i.e. Af = Af^ ® . . . ® Af^ N \ simply 
replace Af„]l % with Af^U in (§§).) 



IV. CONCATENATED CODES 

We now consider concatenated codes [|[ |). We first 
describe the procedure for constructing such codes, and 
then show how the coding maps f2 c make the calculation 
of the effective channels for such codes straightforward. 



A. Constructing Concatenated Codes 

We now describe how two codes may be concate- 
nated to form a larger code; the procedure is depicted 
in Figure |[ Let the two codes be an M-qubit code 
C out = (£ out ,V out ), called the outer code, and an TV- 
qubit code C in = (£ in ,X> in ), called the inner code. A 
logical qubit p is encoded first using the outer code 
C out , yielding the M-qubit state £ out [p ]- Each of these 
qubits is then encoded by the inner code; i.e., the map 
£ in ® . . . ® £ in = (£ in )® M acts on £ out [p }- The compo- 
sition of these encodings forms the encoding map for the 



£ = (£ m ) 



£ c 



(26) 



The M sections of the register encoding each of the M 
qubits in £ out [/9o] are called blocks; each block contains 
N qubits. After the encoding, a noise process Af acts on 
the entire MTV-qubit register. 

A simple error-correction scheme (and one that seems 
reasonable for use in a scalable architecture) coherently 
corrects each of the code blocks based on the inner code, 
and then corrects the entire register based on the outer 
code. I.e., the decoding map for the concatenated code 
is given by 



V = V out o V 



(27) 



We denote the concatenated code (with this correction 
scheme) by C out (C in ) = (£,T>); note that C out (C in ) is 
an MA^-qubit code. 



B. Effective Channels for Concatenated Codes 

Suppose that we have computed the effective channel Q 
due to a code C m = (£ ln ,T> ln ) with some noise dynamics 
Af, and wish to consider the effective channel Q resulting 
from the concatenated code C out (C ln ). We assume that 
each TV-bit block in the register evolves according to the 
noise dynamics Af and no cross-block correlations are in- 
troduced, i.e. that the evolution operator on the MTV-bit 
register is 



Af =Af®Af( 



)Af = Af® M . 



(28) 



By definition, we have Q = V o Af o £ . Substituting (|2 
?]) and (|28|) into this expression yields 



Q = V out o T> imM o Af® M o £™®M £out 
= V out o (V in 0JV0 £ in )«> M o £ out 
= £> out o Q® M o £ out 



(29) 



where we have used Q = T> m oAfo£ m . This result makes 
sense: each of the M blocks of TV bits represents a single 
logical qubit encoded in C ln , and as the block has dy- 
namics Af, this logical qubit 's evolution will be described 
by Q . Comparing with the definition of the coding map 
(|19|), we then have 



g = n c ° u \g) 



(30) 



Thus given the effective channel for a code C m and an 
error model, the coding map fl c makes it straightfor- 
ward to compute the effective channel due to the con- 
catenated code C out (C in ). 

Further, suppose that the original noise model Af 
had the form of uncorrelated noise on single physical 
qubits, as given by <^). Then Q = fl cia (Af w ), and 



£ 



V 



FIG. 2: The error-correction process for the concatenated code C° ut (C in ) = (S,V); here C° ut = (£° ut , V out ) 
code and C" n = (£ ln ,T> ln ) is a five-qubit code. The noise process Af acts on the entire 15-qubit register. 



is a three-qubit 



so Q = £l c ' out (£! c " n (A/'( 1 ))). We may therefore conclude 
that composing coding maps gives the coding map for 
the concatenated code, i.e. 



C" 



(31) 



More generally, we may characterize both the finite 
and asymptotic behavior of any concatenation scheme 
involving the codes {Ck} by computing the maps Vt Ck . 
Then the finite concatenation scheme Ci(C2(. . . C n . . . )) 



is characterized by ^c 1 (c 2 (...c n ...)) _ ^c lo ^c 2( 



We expect the typical fl c to be sufficiently well-behaved 
that standard dynamical systems methods fllPf will yield 
the i — > oo limit of (f2 ) ; one need not compose the 
(Or) explicitly. In section VII, we will consider such 



asymptotic limits in more detail. 



V. DIAGONAL CHANNELS 

As an application of the methods presented above, we 
will consider the commonly-considered error model in 
which each physical register qubit is subjected to the 
symmetric depolarizing channel These single-qubit 
noise dynamics are given by the master equation 

g = \c x \p\ + 1 -C Y [ P \ + \Lz[p\ (32) 

where for any linear qubit operator c the Lindblad deco- 
herence operator C c is given by 

1 t 1 

and 7 is a measure of the noise strength. This master 
equation is easily solved, yielding a qubit channel with 
matrix representation 



C c [p] = epe 1 — -c^cp — -pc^c 



(33) 



M dep = 



1 

e" 7t 








-ft 







-ft 



(34) 



Before calculating effective channels due to this error 
model, it will be useful to discuss the more general set 
of channels whose matrix representation is diagonal. As 
we will see, these channels correspond to single-bit Pauli 
channels, and will allow us to demonstrate the power of 
the techniques developed above. 

Consider a qubit channel given by a diagonal matrix 
J\f^. From trace preservation Afjf 1 = 1, so let the chan- 
nel with Afxx = x i Nyy — V an d ■N'zz - z ^ denoted 
[x, y, z] for compactness. (Thus the depolarizing channel 
@ is given by [e" 7 *, e" 7 *, e" 7 *].) In (llj it is shown 
that complete positivity of such a channel requires 



-x ■ 
x ■ 
x ■ 



y + z < 1 
y + z < l 
y-z < 1 



(35) 



V 



z < 1. 



Now consider the single-bit Pauli channel in which the 
transmitted state is subjected to the Pauli operators X, 
Y, and Z with exclusive probabilities px, Py, and pz, 
i.e. 

P -> (1 - Px ~ PY ~ Pz)p 

+ Px XpX + PY YpY + pzZpZ. (36) 

It is easy to show that this channel has the diagonal ma- 
trix representation 

[l-2{p Y +Pz), l-2(p x +pz), l-2(p x +Py)}, 

(37) 

and so any Pauli channel is a diagonal channel. The 
converse is also true: choosing p x = p y = 

l-x+y-z ^ and ^ _ i-x-y+z yjgi^g tne channel [x,y,z], 
and the complete positivity constraints (|3^) yield the 
standard probability rules px,PY,Pz > and px +Py + 
pz < 1. Thus any diagonal channel may be realized as a 
Pauli channel. Pauli channels are among the most com- 
monly considered error models in the literature, and we 
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will restrict our attention to diagonal channels for the 
remainder of this work. 

The effect of a diagonal channel on a qubit is simple 
to interpret: we have (X) f = x(X)q, (Y) f = y(Y)o, and 
(Z) f = z(Z)q. Thus the X, Y, and Z components of po 
decay independently, and we may therefore speak of the 
deco herence of (X), (Y), and (Z). Recalling from section 
[IB that the fidelity of a pure state p through a qubit 



channel Q is given by ^p Qp, the respective fidelities of 
X-, Y-, and Z-eigenstates through the channel are ^(1 + 
x), and h(l + z). More generally, the fidelity of a 

pure state (requiring (X) 2 + (Y) 2 + (Z) 2 = 1) is given by 
i(l + x{X) 2 +y(Y) 2 + z{Z) 2 ). A common figure of merit 
for a channel is the worst-case fidelity of a pure state, 
which for a diagonal channel is i(l + min(a;, y, z)). Thus 
if for a given error model a code C yields an effective 
channel [x, y, z] and a code C yields an effective channel 
[x' , y', z'], we say that C outperforms C if min(a;, y, z) > 
min(a;', y' ', z'). 

Many commonly considered codes are stabilizer codes 
§, Hi, which are designed to detect and correct Pauli 
errors; it would therefore not be so surprising if the cod- 
ing maps for such codes were particularly well-behaved 
when acting on a Pauli channel. In fact, as proved in Ap- 
pendix W, if C is a stabilizer code and Af^ is diagonal, 
then Sl^A^ 1 ') is also diagonal. Thus just as arbitrary 
codes act as maps on the space of qubit channels, sta- 
bilizer codes act as maps on the space of diagonal qubit 
channels. 



VI. EXACT PERFORMANCE FOR SEVERAL 
CODES OF INTEREST 

We will now present the effective channels for several 
codes of interest under diagonal error models. The codes 
considered here may all be formulated as stabilizer codes; 
thus, as described in the previous section, the effective 
channels will also be diagonal. The diagonal elements 
of the effective channel Q — fl c ([x, y, z]) may be calcu- 
lated either using the coding map methods presented in 
section pi], or using the stabilizer formalism as shown 
in Appendix which may be computationally advan- 
tageous. For each code, we will compute the effective 
channel for a general diagonal error model [x,y,z], and 
then interpret the results for the symmetric depolarizing 
channel Af? cp = [e" 7 *, e" 7 * e" 7 *!. 



The bitflip code first mentioned in section II A is a sta- 
bilizer code; letting f2 bf denote the corresponding coding 
map, we find 



n bt ([x,y,z]) 



[x 3 



x y 



1„3 3. 



(38) 



As the bitflip code is only a three-qubit code, it is not 
unreasonable to check this result with more conventional 
methods, e.g. by counting bitflip and phaseflip errors, 
or by working in the Heisenberg picture to compute the 
evolution of the relevant expectation values. However, 
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FIG. 3: The effective channel [x bt (t),y b( (t), z h{ (t)] due to the 
bitflip code under the symmetric depolarizing channel. The 
respective fidelities of X, Y and Z eigenstates for correction 
performed at time t are given by |(1 + x hi {i)), |(1 + y ht (t)), 
and \(l + z hi {t)). 

for larger codes such computations will rapidly become 
unmanageable . 

To examine the bitflip code acting under the symmet- 
ric depolarizing channel, define [x hl (t), y b! (t), z hi (t)] = 
ft bf (A/" t dop ); the functions x M , y hf , and z M are plotted in 
Figure [| along with e -7 * (describing the decoherence of 
the physical qubits) for comparison. We see that z hi (t) > 
e -7 *, and thus the decoherence of (Z) is suppressed by 
the bitflip code. However, x hf (t) = y hl (t) < e~ 7 *, and 
thus the decoherences of (X) and (Y) are increased by 
the bitflip code. 



More generally, for any < x, z < 1 we have ■ 



z and x 3 < x, so for any physical channel in this regime 
the bitflip code always suppresses decoherence of (Z) and 
increases decoherence of (X) . Decoherence of (Y) is sup- 
pressed when x > ^/| and < y < y/3x 2 — 2, and in- 
creased for all other positive values of x and y. We may 
therefore conclude that under a general Pauli channel 
the bitflip code increases the fidelity of some transmitted 
states at the expense of others, and thus the bitflip code 
is outperformed by storing the logical qubit in a single 
physical bit. 

However, as the bitflip code is designed to only pro- 
tect against physical bitflip (X) errors, it should not be 
expected to perform well in the presence of arbitrary 
Pauli errors. If we consider physical channels with only 
X errors, we find that the bitflip code suppresses de- 
coherence of all encoded states. More precisely, sup- 
pose that the physical qubits are evolving via a Pauli 
channel ( ^6|) with only X errors, i.e. py = Pz = 0. 
Then [x,y,z] = [1,1-2^,1-2^], and Q hi ([x, y, z]) = 
[1,1— \p\ + \p\-, 1 — \p\ + \Px\- Thus we have repro- 
duced the usual result of a leading-order analysis: the 
bitflip code suppresses decoherence due to X errors to 
order p 2 x . 
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Now consider the three-qubit phaseflip code M, with 
encoding |±) i-> | ± ±±) for |±) = -^(|0) + |1)). This 

code is completely analogous to the bitflip code, detecting 
and correcting single phaseflip (Z) errors instead of single 
bitflip (X) errors. The phaseflip code's coding map Q pf 
is exactly the same as that of the bitflip code, with the 
role of X and Z interchanged: 



n* l ([x,y,z])= [|x-ix 3 , \z\-\y\ z 3 



(39) 



The concatenation phaseflip(bitflip) yields the Shor nine- 
bit code |], | with encoding |±) i-> ^(|000) ± |111))® 3 . 

Thus Q sho1 = Q pi o fl hf . Evaluating this composition 
using the coding maps (p8|) and (p^ 



Q^([x,y,z}) 



n pi (n hi ([x, y ,z])) 

[P(x),Q(x,y,z),R(z)] 



(40) 



where 



P(x) 
Q(x,y,z) 

R(z) 



3-3 _ 1 T 9 
2 X 2^ 



_l(3^ y _l y3) 3 



(41) 



(The combinatoric analysis required to reproduce this re- 
sult by counting bitflip and phaseflip errors would be 
quite tedious!) 

To examine the Shor code acting on the symmetric 

z s ~ 
and 



depolarizing channel, let [:r shor (i) 



y Shor (t), 



„Shor „.Shor 



v Shor 



arc 



n Shor (A/7 ep ); the functions ar— , y" 
plotted in Figure We see that for short times (or equiv- 
alently, weak noise-strength 7) , the Shor code suppresses 
decoherence of (X), (Y), and (Z). For long times, how- 
ever, the code increases the decoherence of all three ex- 
pectation values, and as z Shm (t) > x shol {t) > y shm (t), in 
an intermediate regime the code suppresses the decoher- 
ence of some of the expectation values while increasing 
that of others. Thus to suppress the decoherence of an 
arbitrary logical state, correction needs to be performed 
at a time t when y shol (t) > e -7 *. 

Above we defined the phaseflip code by the encoding 
|±) 1— > I ± ±±); we could have also used the encoding 

given by |0) h-> | + ++>, |1> h-> | ). Call the code with 

this encoding phaseflip', with coding map O pf . As this 
modification of the phaseflip code simply interchanges 
the encoded X- and Z-eigenstates, the new effective 
channel is simply that of the original phaseflip code with 
the effects of the channel on the X and Z components 

x,y,z\) = [z 3 ,|z 2 ?/-i?/ 3 ,fa;- 
Wc could then use this version of 



of po interchanged: S! pf ([. 



(compare to (|39|)). 



the phaseflip code to define an alternative version of the 
Shor code with the encoding |0) i-> ^=(|000) + llll))' 83 , 

|1) 1 > -1=(|000> - |111))® 3 . Call this code Shor', with 
corresponding coding map 17 shor = Vt v o f2. We find 
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FIG. 4: The effective channel [x ShoT (t) , y ShoT (t) , z ShoT (t)] due 
to the Shor code under the symmetric depolarizing channel. 

with the polynomials P, Q and R defined by (|4l]). Com- 
paring to (40), we see that again this modification of the 
Shor code simply interchanges the effect of the channel 
on X and Z components of po. Assuming that the en- 
coded logical states are randomly distributed (as opposed 
to always sending Z-eigenstates, for example), the choice 
of using the Shor code or the Shor' code is simply one of 
aesthetics: the effective channels are identical up to in- 
terchange of the decoherence of (X) and (Z). However, 
as we will see in the next section, this choice does have 
an impact when these codes are concatenated. 

For comparison, we consider two other stabilizer codes 
of interest. The Steane code [fl f| is a seven-bit code de- 
signed to correct errors consisting either of a Pauli error 
(X, Y, or Z) on a single qubit of the codeword, or of an 
X and a Z error on separate qubits. We find 



Stcanc 




with 

S(x) 
T{x,y,z) 



{[x,V, z}) = [S(x),T(x,y,z),S(z)} 



(43) 



7 3 _ 3 T 7 

y + & 7 



16 



(x 4 + z 4 )y 3 + fx 2 yz 2 . 



(44) 



n Sieaae (N? ep 

find that the functions x stcanc , y Htoane and z st =<"" arc 



Let [x stcanc (t), y stcanc (i), z Stoanc (i)] 

„Stcanc „ .Stcanc 



t >, we 

Stcanc 



qualitatively similar to the analogous functions of the 
Shor code. If they were plotted in Figure |], these func- 
tions would be interspersed between the plotted curves: 
for all values of t > 0, we have z Shor > z Stcanc = x stcanc > 
^Shor > y steane > y Shor Though the Shor code more ef- 
fectively suppresses decoherence for logical Z-eigenstates, 
the Steane code performs better in the worst case (Y- 
eigenstates), and thus outperforms the Shor code. 

The Five-Bit code [§ |l|, corrects Pauli errors on 
a single qubit of the codeword. We find 



Shor' 



([x,y,z]) = [R(z),Q(x,y,z),P(x)}. 



(42) 



n Five ([x, y, z]) = [U(x, y, z), U(y, z, x), U(z, x, y)} (45) 
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with 



U(x,y,z) = \x{y 2 + z 2 ) 



\xy 2 z 2 



1-5 
4 x . 



(46) 



Five/- 



Letting [x 

„Five „,Five 



||,ivo (t),z Fivc (i)] = fi Five (AA t dcp ) yields 
as expected from the symmetries of 



(t),y 

x — = r — = z Five 

the code and of the map fi Flve . Thus the fidelity of a state 
through this channel is independent of the state. These 
functions also have qualitatively similar behavior to those 
plotted in Figure f| and for t > we have z ShoT > z Fivc > 
z stoano > ^Shcff Thug ^ p ive _ Bit code outperforms 

both the Shor and Steane codes. 



VII. EXACT PERFORMANCE AND 
THRESHOLDS FOR CERTAIN 
CONCATENATION SCHEMES 

We now consider the effective channel due to fami- 
lies of concatenated codes under the symmetric depo- 
larizing channel. First, consider the Shor code concate- 
nated with itself £ times, denoted by Shor*". From section 
IV, we know that the coding map for this code is given 



by fi( Shor ) = tt shm o . . . o O shor = (n shm ) e . As fi Shor 
takes diagonal channels to diagonal channels , the effec- 
tive channel due to Shor is also diagonal. Let 



{x e (t),y e (t),z e (t)} = (n» h ° r yWt P ), 



(47) 



)Shor 



which may be calculated using the polynomials of 
given in (^Tj). 

The functions zt(t) for < £ < 4 are plotted in 
Figure |[ We immediately observe that as £ — ► oo 
the functions zt(t) approach a step function. Denoting 
the step function's time of discontinuity by t* Zl we have 
zt(t) -> 9{t* z - t) where 



9{x) = 



x < 

1 x > 



(48) 



For t <t z , each layer of concatenation decreases the (Z) 
decoherence, yielding perfect preservation of the encoded 
(Z) information in the infinite concatenation limit. How- 
ever, for t > t* z , the (Z) decoherence increases. Thus in 
the infinite concatenation limit, the code will perfectly 
protect (Z) of the logical qubit if correction is performed 
prior to t* z ; if correction is performed after this time, all 
(Z) information is lost. 

Similarly, the functions X((t) and ye(t) approach step 
function limits as £ — > oo; call the discontinuous times 
of these step functions t* x and t Y . To perfectly protect 
an arbitrary state in the infinite concatenation limit, cor- 
rection must be performed prior to ith = min(t^ , t Y , t* z ). 
We call tth the storage threshold. (We use the term "stor- 
age threshold" to indicate that the threshold takes into 
account only noise in the register, rather than gate or 
measurement errors also considered in fault-tolerant set- 
tings.) We now show how the coding map Q shm may be 
used to find this threshold. 
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FIG. 5: The functions zt, where [xe(t), ye(i), ze(t)] is the ef- 
fective channel for £ concatenations of the Shor code under 
the symmetric depolarizing channel. 




FIG. 6: The function R(z) (plotted as the thick curve). Ob- 
serve that the map z \— > R(z) has fixed points at 0, 1, and z*. 
The arrows depict the iteration of this map pushing points in 
the interval (0, z*) toward and points in the interval (z* , 1) 
toward 1. 

Observe in Figure || that the plots of Zt(t) all intersect 
at a point (jt* z ,z*). Writing f2 Shor in the form (Eoh, we 
have ze +1 (t) = R(ze(t)). The function R(z) is plotted 
in Figure ^. We see that the map z i— » R(z) has fixed 
points at 0, 1, and a point z* » 0.7297. (We find z* 
by numerically solving z = R{z) on the interval (0,1).) 
Iterating the map pushes points in the interval (0, z*) 
toward 0, and pushes points in the interval (z* , 1) toward 
1 . In the language of dynamical systems |l(J , and 1 are 
attracting fixed points, and z* is a repelling fixed point. 
This behavior leads to the shape of the plots in Figure 



We then invert e 



-7*z 



to find -yt* z = 0.3151. 



The function P(x) has the same qualitative behavior on 
(0, 1) as R(z), so we may similarly find x* w 0.9003 and 
-it\ fa 0.1050. 
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Code 


Shor 


Shor' 


Steane 


Five-Bit 


a 


X,Y 


Z 


X,Y 


Z 


X,Y,Z 


X,Y,Z 




0.1050 


0.3151 


0.1618 


0.2150 


0.1383 


0.2027 


Pth 


0.0748 


0.1121 


0.0969 


0.1376 



TABLE I: Code storage threshold results. (See text for dis- 
cussion.) 

We cannot use the same method to find t Y , as yg + i(t) is 
a function of xi{t), yi(t) and zt(t), not just of yi(t) alone. 
(This problem is evident from plots of the functions yi(t): 
though these functions approach a step function in the 
£ — > oo limit, they do not all intersect at a point as the 
plots of Zi{t) do.) However, we now argue that finding t* x 

xe(t) — > 1 and 



and t* z is sufficient to find t Y 



For t < t* x , 



zt(t) — > 1. Using the complete positivity constraints (35), 
we find that if [x,y,z] is a channel, x = z — 1 implies 
y = 1. Since the space of channels [x, y, z] is closed and 
bounded (it consists of the boundary and interior of a 
tetrahedron in R 3 ), xi{t) — > 1 and Ztit) — > 1 implies 
y^(t) — > 1. Now for t* x < t < tig, x e (t) — > and z e (t) — > 1. 
Using the complete positivity constraints (p5[), we find 
that if [x, y, z] is a channel, x — and z = 1 implies y = 0. 
Thus we may conclude that for these values of t, ye{t) — > 
0. We now have y^(i) — > 1 for t < t* x , and yt(t) — ► 
for t > t^r, thus we conclude t Y — t* x . More generally, 
if we know t 



x 



and t* z , then fcy- 



is given by min(t x ,t z ). 
We may therefore conclude that jt Y as 0.1050, and so 
jt t h as 0.1050. (The value of t Y could also be obtained 
from the dynamics of the polynomial maps P, Q and 
R without making reference to the complete positivity 
constraint, but the method presented here requires less 
argument at ion . ) 

We may also phrase these thresholds in the language 
of finitely probable errors. Consider the symmetric Pauli 
channel given by ( |36| ) with p x = py = pz = §■ 
This channel subjects a qubit to a random Pauli er- 
ror with probability p, and is described by Af Faull (p) = 
[1 — |p, 1 — |p, 1 — |p] . Observe that the symmetric Pauli 
channel and the symmetric depolarizing channel are re- 



dep 



Thus in the infinite 



lated by A/" Pauli ( § ( 1 - e -T ) ) 
concatenation limit with AA Paull (p) acting on each regis- 
ter qubit, the (er) of the logical qubit will be perfectly 
protected if p < p* — |(1 — e~ 7 *<^). Define the threshold 
probability p th = mm{p x ,p Y ,p z }; for p < p th , all en- 
coded qubits are perfectly protected in the infinite con- 
catenation limit. Values for 7^* and pth appear in Table 

s 

We now use similar methods to derive thresholds for 
the Shor', Steane, and Five-Bit codes presented in the 
previous section. First, consider the Shor' code. Let 
[x' e (t),y' e (t),z' e (t)} = (r» ShOT ') £ (<° P )- The y' e (t) approach 
a step function as I — > oo, but x' t {t) and z' e (t) approach 
a limit cycle of period 2: we find that x' 2i and z' 2i+1 both 
approach 0{t\ — t) for some value of t*, while x' 2l+1 and 
approach 0(t 2 — t) for some distinct value of t 2 . From 



-21 



the form of Vr given in ([42j), we see that x' e+1 (t) is 
a function of z' e (t), and z' g+l (t) a function of x^it), so 



it is not so surprising that the sequence z' Q , x\, z' 2 , x' 3 , . 
converges and the sequence x' , z[ , 



converges. 

To find the threshold, we simply consider the sequence 



b 2' *3) ' 



of channels 



[x 2e , y 2£ , 



z' 2l \, generated by the map (f2 Shor ) 5 



From (42) we see that x^+i) = R(P( x 2t)) an( i z 2(i+i) ~ 
P(R(z' 2e )). Thus to find the values t* x , t Y and tig, we 
find the fixed points of the maps x > R(P(x)) and 
z i ^ P(R(z)), and proceed as with the Shor code. As 
shown in Table |, we find that, compared to the Shor 
code, the Shor' code has greater values for t* x and t Y , 
and a lesser value for t* z . As the threshold tth is given 
by the minimum of these three values, the Shor' code 
outperforms the Shor code in the infinite concatenation 
limit. 

The map f] stcanc ; given by ([i"3|), has the same form 
as the Shor code map (|40|), and therefore we can use 
the same methods to find the Steane code thresholds. 

has a different form. 



The map fl c , given by 
However, as Af t cp has the symmetric form [ 1 and 
fl e preserves this symmetry by taking [ to 
[U (x, x,x),U [x, x,x),U [x, x, x)], we may find ti x = t Y = 
ti z simply by finding the fixed point of x i— ► U(x,x,x). 
Results are summarized in Table [jj. We find that the 
Five-Bit code has the largest threshold, and therefore 
the best performance in the infinite concatenation limit. 
It is interesting to note that the Shor' code outperforms 
the Steane code in the infinite concatenation limit, even 
though the opposite is true for only one layer of each 
code. 

We conclude our discussion of the thresholds by com- 
paring the exact values of pth to those calculated with 
traditional leading-order techniques (e.g. in |3|). First 
consider the Five-Bit code. Under the symmetric Pauli 
channel A/" Pauh (p), a physical qubit is unmodified with 
probability 1 — p. The Five-Bit code perfectly protects 
its encoded information if no more than one of the five 
physical qubits are subjected to a Pauli error. Under 
A/" Pauh (p) , the probability of no errors on any physical 
qubit is (1 — p) 5 , and the probability of exactly one error 
is 5p(l — p) A . We then assume that all greater-weight er- 
rors are uncorrectable, and find that the probability of a 
correctable error is (1— p) 5 + 5p(l— p) A = 1 — I0p 2 +O(p 3 ). 
The threshold value pth is the value of p at which the sin- 
gle physical qubit and the encoded information have the 
same probability of error. Thus to estimate the thresh- 
old we solve 1 — 10p 2 = 1 — p, yielding p t h — jq- Thus 
the leading-order calculation underestimates the actual 
threshold (0.1376) by 27%. (The assumption that all 
errors of greater weight are uncorrectable assures that 
the approximation underestimates, rather than overes- 
timates, the threshold.) The Steane code corrects all 
weight-one errors, and weight-two errors consisting of an 
X on one bit and a Z on another bit. A similar calcu- 
lation finds the probability of a correctable error to be 
1 - f p 2 + 0(p 3 ), yielding p th m 0.0612, a 37% under- 
estimate. The Shor code corrects all weight-one errors, 
and weight-two errors such that any X and Y operators 
occur in different blocks, and any Y and Z errors occur 
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in the same block. The probability of a correctable error 
is found to be 1 - 16p 2 + 0(p 3 ), yielding p th = 0.0625, 
a 16% underestimate. The analysis is exactly the same 
for the Shor' code, yet the Shor' code has a very different 
threshold; in this case, the leading-order result underes- 
timates the threshold by 44%. 



VIII. CONCLUSION 

We have shown how a code's performance for a given 
error model can be described by the effective channel 
for the encoded information. The methods presented 
for calculating the effective channel have allowed us to 
find the performance of several codes of interest under 
single-bit Pauli channels, and further have allowed us to 
find thresholds describing these codes' asymptotic lim- 
its of concatenation under the symmetric depolarizing 
channel. Though we chose to restrict our attention to 
diagonal channels, these methods can be applied to any 
uncorrelated error model (e.g. the amplitude-damping 
channel ||, which is non-diagonal), and will substantially 
simplify the exact analysis of code performance in these 
more general settings. 

We believe that this effective channel description of 
code performance may be useful in other contexts as well. 
For example, this work could be extended to take ac- 
count of encoding and decoding circuit errors, thereby 
providing a method for calculating exact fault-tolerant 
thresholds. Also, by providing a comprehensive method 
for describing the performance of a quantum code with- 
out reference to a particular error model (e.g. bitflip and 
phaseflip errors) perhaps these methods will allow us to 
address open questions such as the optimal code for a 
given error model, and the quantum channel capacity. 
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APPENDIX A: STABILIZER CODES AND 
DIAGONAL CHANNELS 

In this appendix we consider the effective channel Q = 
H, c (Af^) when AA 1 ) is diagonal and C is a stabilizer 
code. We show that Q is also diagonal, and show how the 
stabilizer formalism facilitates its calculation. The reader 
unfamiliar with stabilizer codes is directed to || for an 
introduction, and to (l2) for a more complete discussion. 

Since AA 1 ) is diagonal, the terms Nvljn in the expres- 
sion for the effective channel (p5|) vanish for ^ /i.;. 



Thus we have 



JV 



(Al) 



dramatically simplifying the calculation of Q. The coef- 
ficients a", i and /3f x are defined in terms of the E„i 



and D a operators in (]2(]) and (|TJ); to calculate these 
operators we now consider the code C in more detail. 

Let C be a stabilizer code given by stabilizer S = 
{S k } C ±{7, X, Y, Z}® N , storing one qubit in an iV-qubit 
register. The stabilizer S defines the codespace, and 
the logical operators I,X,Y,Z C ±{I,X,Y,Z}® N de- 
termine the particular basis of codewords |0), |T). Recall 
that the E a > operators act as \a' on the codespace and 
vanish elsewhere. It can be shown that Pq — jh^ Yl k St 
acts as a projector onto the codespace, and by definition 
the logical operators a 1 act as a 1 on the codespace. Thus 



E„ 



(A2) 



will act as ha' on the codespace and vanish elsewhere. 

As an example, consider the bitflip code introduced in 
section [I A. The bitflip code may be specified as a sta- 



bilizer code, with S = {III_,ZZI,IZZ,ZIZ}, I = III, 
X = XXX, Y = —YYY, Z = ZZZ. The above expres- 
sion reproduces the expressions for the E a i presented in 
([22l). Without the stabilizer formalism, deriving ( ^2|) is 
an exercise in expanding projectors in a basis of Pauli op- 
erators; with this method the computation is very simple. 

We now construct the D a operators for the stabilizer 
code. As in section [I A, let {Pj} be the projectors de- 
scribing the syndrome measurement. For a stabilizer 
code, the recovery operators Rj are each chosen to be 
a Pauli operator taking the space projected by Pj back 
to the codespace. Consider the expression for D a given 
by (0); substituting in the expression ( |A2| ) for E a , we 
have 



(A3) 



Now because Rj, Sk and a are all Pauli operators, they 
either commute or anti-commute. For two Pauli oper- 
ators V and W, let T)(V, W) = ±1 for VW = ±WV. 
Commuting the Rj to the left in the above expression 

and noting that RjRj = 1, 



D <y = jk\ Y^k.j -ri{S k ,Rj)ri(Rj,a)S k cT 



(A4) 



with fka = J2j r l(^k,Rj)v{RjT a )- Again, as an example 
consider the stabilizer definition of the bitflip code. The 
recovery operators are III, XII, IXI, and IXI. Eval- 
uating the above expression for D a yields the previous 
result of (|23|). 
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Using the expressions (A2) and (A4) for the E a i and 
D a operators in the stabilizer formalism, we we will now 
find the coefficients otf and Pf Vi \ as defined in (gC ) 
and (2lh. Since SV-cr is a Pauli operator, the sums (|A2) 
and (A.4) are expansions of these operators E a i and D a 
in the Pauli basis; if we were to write down these sums 
explicitly for a given stabilizer code, the coefficients a 
and 8 could be read off immediately, e.g. from ( p2|) and 

This approach may be formalized as follows. First, 
note that Sk and cr are both Hermitian Pauli operators, 
and they commute; therefore their product is also a Her- 
mitian Pauli operator, i.e. S k a € ±{I, X, Y, Z}® N . For 
any operator V = ±/xj ® . . . ® /x^r with yn E {I, X, Y, Z}, 
let (/>(V) =/Xi®...®/ijv, and let a(V) € {0, 1} such that 
7 = (-l]< y V(7). 
re-write (|A2j) and (t 



Then, using jS 1 ! = 2 W , we may 



-Do- = 



k 



^_^a(S'fcO-')_ 



hr<j>{S k CF') 



k 



{-l) a{s * 9) ±h^{S k o). 



(A5) 



(A6) 



Comparing (A5) with the definition of 



in (20) 



we see that each term of the sum over k contributes to 
a single coefficient o^/ Sk9 ,-., as 7pr<fi(Sf.a') is of the form 

® ■ ■ ■ <8> 4a*jv)- Similarly, each term in (A6) con- 
tributes to a single coefficient P^g kS y Lemma 2 of Ap- 
pendix [B| shows that (fi(Skcr) ^ 4>(S)~' a) unless k = k' and 
a = a' . Thus each term in (|A5|) contributes to a distinct 



coefficient 



l 0(Sk<7') 



and each term in (|A6|) contributes to 



a distinct coefficient /3T 



0(S fc <r)' 



We may therefore simply 



read off the coefficients from (A5) and (A6), yielding 



(Sua' 



a(S k a') 



PI 



<P{S k 9) 



(-1) 



= (-i) a(Sfcff) ^JW, 



(A7) 
(A8) 



and all other ot a s -, and 87 , vanishing. 

We now evaluate (J = fl c (Af^ 1 ') where A/"*- 1 ^ = [x, y,z] 
using (Al). The only non- vanishing terms Pf^.x occur 
when fii ® • • ■ <8 A*JV = for some fc and cr, and 

the only the only non- vanishing terms a^.j occur when 
[i\ ® . . . ® jiijv = 4>(Sk&') for some fc and cr'. Thus the 
coefficients /?? 
^at = <j)(S k d) 



of ( Al) will vanish unless |«i<8>. . .0 
4>{Ski(j') for some fc and fc'. As proved 
in Lemma 2 of Appendix |B], this cannot happen when 
a ^= a' . Thus all the matrix elements Q a <j' vanish when 
a 7^ cr', i.e. (J is diagonal. 

Having demonstrated that the coding map Q of a sta- 
bilizer code C takes diagonal channels to diagonal chan- 
nels, and because Qu = 1 from trace preservat ion, we 
need only compute Qxx, Qyy, and Qzz using (Al) to 
find Q = fl c ([x, y, z]). These computations can be per- 
formed using the methods of section [II , but we conclude 
this section by expressing these elements using the stabi- 
lizer formalism, which may be computationally advanta- 
geous. 



Consider the diagonal terms Q aa given by ( |Al[ ). We 
need only sum over the non-vanishing coefficients a and 
8, which are given by (A7) and flA8|). Substituting in 
these expressions yields 



(A9) 



where 4>i{V) denotes \ii for (f>(V) — fii (8 • • • <8 ^N- Now 
as = [x,y,z], the product of the matrix elements 

of AA 1 ) in the previous expression is simply a product 
of x's, y's and z's; each factor appears as many times 
as (respectively) X, Y and Z appear in (f)(S k a). Letting 
w a (p) denote the a- weight of a Pauli operator p, e.g. 
wx(XYX) = 2, we have 



x w x {Sk<7) yWY {Sk<7) z wz{S k cr) 



(A10) 



APPENDIX B 



This appendix contains lemmas deferred from previous 
sections. 

Lemma 1: The decoding operation V given by (^) is 
a quantum operation. 
Proof: From (Q) we have 



(Bl) 



To prove that I? is a quantum operation, we must show 



(B2) 



where 1 is the identity on the register space. As we as- 
sumed that 1Z maps all states into the codespace, we 
can choose the operators Aj to only have range on the 
codespace. With such a choice, A'.Aj = Aj(|0)(0| + 



|r)(l|)A/ = A]BB^Aj. We therefore have 
5>}B)(B%-) = 5>^- 



(B3) 



From (|3|) we have J^j AjAj — 1, and so I? is a quantum 
operation. ■ 

Lemma 2: For a stabilizer code with stabilizer {SjA 
and logical operators {cr}, and (j> defined in Appendix 0, 
4>(Skd z ) ^ (/i(S'fc'cr') unless k = k' and cr = cr'. 

Proof: Suppose we have c6(Sfccr) = ci(S'fc'cr'); then 
Sfccr = ztSk'd- 1 . As the stabilizers Sk and Sk' act triv- 
ially on the codespace, Sfccr and ±Sk'Cf' act respectively 
as cr and ±cr' on the codespace. Thus we must have 
cr = ±cr', which requires a = a' and the ± sign be pos- 
itive. We now have Sk& = Sfc'cr; right-multiplying by cr 
yields Sk — Sk> , and thus k = kl . ■ 
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